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We offer a new, physically transparent argument for the existence of the critical, universal 
maximum matter density in loop quantum cosmology for the case of a flat Friedmann-Lemaitre- 
Robertson- Walker cosmology with scalar matter. The argument is based on the existence of a sharp 
exponential ultraviolet cutoff in momentum space on the eigenfunctions of the quantum cosmological 
dynamical evolution operator (the gravitational part of the Hamiltonian constraint), attributable 
to the fundamental discreteness of spatial volume in loop quantum cosmology. The existence of the 
cutoff is proved directly from recently found exact solutions for the eigenfunctions for this model. 
As a consequence, the operators corresponding to the momentum of the scalar field and the spatial 
volume approximately commute. The ultraviolet cutoff then implies that the scalar momentum, 
though not a bounded operator, is in effect bounded on subspaces of constant volume, leading to 
the upper bound on the expectation value of the matter density. The maximum matter density is 
independent of the quantum state essentially because of the linear scaling of the cutoff with volume. 
These heuristic arguments are supplemented by a new proof in the volume representation of the 
existence of the maximum matter density. The techniques employed to demonstrate the existence 
of the cutoff also allow us to extract the large-volume limit of the exact eigenfunctions, confirming 
earlier numerical and analytical work showing the eigenfunctions approach superpositions of the 
eigenfunctions of the Wheeler-DeWitt quantization of the same model. We argue that generic (not 
just semiclassical) quantum states approach symmetric superpositions of expanding and contracting 
universes. 
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I. INTRODUCTION 

Loop quantized cosmological models generically predict that the "big bang" of classical general relativity is replaced 
by a quantum "bounce" in the deep-Planckian regime, at which the density of matter is bounded by a maximum 
density, typically called the "critical density" Pcrit- (See Ref. [TJ[2] for recent reviews of loop quantum cosmology [LQC] 
and what is currently known about these bounds in various models, as well as references to the earlier literature.) 
In most models, the value of this critical density is inferred from numerical simulations of quasiclassical states. 
So far it has been possible in only a single model - the exactly solvable loop quantization, dubbed "sLQC" 
of a flat Friedmann-Lemaitre-Robertson- Walker cosmology sourced by a massless, minimally coupled scalar field - 
to demonstrate analytically the existence of pcrit for generic quantum states. In this model, it was shown that 
Pcrit ~ 0.41pp, where pp is the Planck densityj^ As in other models, the bound for the density was found first 
numerically in Refs. [SIS]. This value was then confirmed and given a clean analytic proof in Ref. [3 . 

In this paper we offer a new demonstration of the existence of a critical density in this model with the hope of 
enriching the understanding of existing results. The argument is rooted in a study of the behavior of the dynamical 
eigenfunctions of the model's evolution operator, the gravitational part of the Hamiltonian constraint, based on an 
explicit analytical solution for these eigenfunctions found recently in Refs. [71[H]. We will show from this solution that 
the eigenfunctions exhibit an exponential cutoff in momentum space that is proportional to the spatial volume. This 
ultraviolet cutoff may be understood as a consequence of the fundamental discreteness of spatial volume exhibited 
by these models. As a consequence of the cutoff, the quantum operators corresponding to the scalar momentum and 
spatial volume approximately commute. The ultraviolet cutoff then implies that the scalar momentum - even though 
its spectrum is not bounded - is in effect bounded on subspaces of constant volume. The proportionality of the cutoff 
to the spatial volume then leads to the existence of a critical density that is independent of the quantum state. 

It has long been understood in the loop quantum cosmology community that the behavior of the eigenfunctions of 
the gravitational Hamiltonian constraint operator is the key to understanding the physics of loop quantum models. 
In particular, the quantum "repulsion" generated by quantum geometry at small volume - leading to the signature 
quantum bounce - was clearly recognized in the decay of the eigenfunctions at small volume in many examples |5l [6] . 
(See also e.g. [5HTT]. among many others.) It was also recognized numerically that the onset of this decay as a function 
of volume depended linearly on the constraint eigenvalue. (See e.g. [HI HH]-) What is new in this work is the shift in 
focus to the behavior of the eigenfunctions as functions of the continuous variable k labeling the constraint/momentum 
eigenvalues. This allows certain insights that may not be as evident considered as functions of the discrete volume 
variable v. From the exact solutions for the eigenfunctions of sLQC, we are able to show a genuinely exponential 
cutoff in the eigenfunctions as functions of k that sets in at a value of k that is proportional to the spatial volume, thus 
confirming and grounding the numerical observations analytically. This, of course, is the same cutoff that manifests as 
the decay of the eigenfunctions at small volume, considered as functions of the volume ~ this is clearly evident in, for 
example. Fig. [3]- but seen from a complementary perspective that is in some ways cleaner because of the continuous 
nature of the variable k. From there, we go on to show how the linear scaling of the cutoff gives rise to the universal 
upper limit on the matter density. To our knowledge, this connection between the linear scaling of the cutoff on the 
eigenfunctions and the existence and value of the universal critical density has not previously been noted. 

Though appealingly intuitive, this argument is essentially heuristic, so we supplement it with a new proof in the 
volume representation of the existence of a Pcrit in this model. (The proof of Ref. [5] is in a different representation 
of the physical operators.) 

Thus we are able to offer a clear physical and mathematical account of the origin and value of the critical density, 
grounded analytically in the exact solutions for this model, that complements and confirms extensive numerical and 
analytic results extant in the literature. This perspective may be of some use in numerical and analytical investigations 
into the existence of a critical density in more complex models for which full analytical solutions are not available. 
We expand on this point in the discussion at the end, after we have developed the necessary details. 

As a by-product of the methods employed to reveal the ultraviolet cutoff on the dynamical eigenfunctions, the semi- 
classical (large volume) limit of the eigenfunctions is also obtained from the exact eigenfunctions. The result confirms 
the essence of the result obtained on the basis of analytical and numerical considerations in Refs. (5l[6j, that the exact 
eigenfunctions approach a linear combination of the eigenfunctions for the Wheeler-DeWitt quantization of the same 
physical model. (See also Ref. [TT], in which a careful analysis of the asymptotic limit of solutions to the gravitational 
constraint arrived at the same result as demonstrated here from the explicit solutions for the eigenfunctions.) The 
domain of applicability of this approximation is described. This result is then used to argue that, in the limit of large 
spatial volume, generic states in LQC - not just quasiclassical ones - become symmetric superpositions (in a precise 



^ The difference of a factor of 1/2 in the value of PcHt quoted here to that in the earhest papers is attributable to the realization in Ref. 
[4] that the "area gap" A of loop quantum gravity should contain an additional factor of 2 for these models; see footnote 1 of Ref. [3]. 
See also footnote Isl 
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sense to be specified) of expanding and contracting universes. The symmetry exhibited in numerical evolutions of 
semiclassical states - see e.g. Refs. [U |3l [3 [6] - is therefore not an artifact of semiclassicality, but a generic property 
of all states in loop quantum cosmology. (Compare Refs. [TTl [T2] for analytic results bounding dispersions of states, 
showing they remain small on both sides of the bounce.) 

The plan of the paper is as follo ws. In Sec.[ll]we summarize the loop quantization of a flat FLRW spacetime sourced 
by a massless scalar field. Sec. 



Ill 



studies the dynamical eigenfunctions e^*^' {i>) of the model in detail, exhibiting various 



explicit forms for the solutions, and works out the asymptotic behavior of the 
A|fc| ^ (The ultraviolet cutoff on the e\f\v) emerges from this analysis in Sec 



in the limitsM S> A|A;| and 



III A 3 



) In Sec. 



IIIB 



the cutoff is 



employed to place bounds on the matrix elements of the physical operators and argue that the scalar momentum is 
approximately diagonal in the volume representation. Section IV applies these results to show that generic states in 
sLQC are symmetric superpositions of expanding and contracting Wheeler-DeWitt universes at large volume. Finally, 
Sec. [V| offers an intuitive argument for the existence of a critical density in this model based on the UV cutoff for the 
eigenfunctions, as well as a new analytic proof in the volume representation. Section [Vl| closes with some discussion. 



II. FLAT SCALAR FRW AND ITS LOOP QUANTIZATION 



In this section we briefly describe the loop quantization of a flat {k = 0) Fricdmann-Robertson- Walker universe 
with a massless, minimally coupled scalar field as a matter source. The model is worked out in detail in Refs. [3l El |6] 
(see also Ref. il3|); see Ref. |7j for a summary with a useful perspective and Refs. [Hll] for recent general reviews of 
results concerning loop quantizations of cosmological models. 



A. Classical homogeneous and isotropic models 

The starting point is a flat, flducial metric qab on a spatial manifold S in terms of which the physical 3- metric is 
given by Qab — o? Qab, where a is the scale factor. The full metric is given by 

gab = -naUb + Qab, (2.1) 

where the normal Ua = —Ndta to the fixed (CtQab ~ 0) spatial slices is given in terms of a global time t and lapse 
N{t), so that a — a{t). 

For the Hamiltonian formulation of the quantum theory spatial integrals over a finite volume are required. We may 
therefore either choose E to have topology T^^ with volume V with respect to Qab, or topology and choose a fixed 
fiducial cell V, also with volume V with respect to qab- The choice plays no role in the sequel and we will proceed in 
the language of the latter choicej^ The physical volume of V is therefore V = a"^ V. 

For a massless, minimally coupled scalar field, after the integration over the spatial cell V has been carried out the 
classical action is 

S^V fdtl-^^-^+'-a^^-]. (2.2) 
J I SttG iV 2 iV J ^ ' 

The classical Hamiltonian is thus 

1 r 2t:GN 21^2! ON 

where Pa and are the canonical momenta conjugate to the scale factor and scalar field. 

Solving Hamilton's equations yields the classical dynamical trajectories, for which is a constant of the motion, 
and 



± , In 
Vl27rG 



(2.4) 



where Vo and c/fo are constants of integration. Regarding the value of the scalar field (j) as an emergent internal physical 
"clock", the classical trajectories correspond to disjoint expanding (+) and contracting (— ) branches. The expanding 



^ For some discussion of this point see Sec. II. A. 1 of Ref. [T]. 
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FIG. 1. Two classical trajectories (Eq. (2.4|) for a massless scalar field in a flat homogeneous isotropic universe are shown. The 



solid (red) curve corresponds to an expanding branch and the dashed (blue) curve to the corresponding disjoint contracting 
branch. The branches are singular in the "past" and "future" given by the internal time 0, respectively. (Figure taken from 
Ref. M-) 



branch has a past singularity (the big bang) in the limit (f> — > — oo, and the contracting branch a future singularity 
(big crunch) as — +oo. (See Fig. [I]) Note that all classical solutions of this model are singular in one of these 
limits. 

Finally, we observe that the matter density on the spatial slices S at scalar field value (f> is given in the classical 
theory by the ratio of the energy in the scalar field to the volume at that (p, p — TabU°'u'', where — [d/drY and 
dT = N dt: 



P\4> 



2V\ 



(2.5) 



B. Loop quantization 

In the quantum theory, following Ref. ^ we will discuss volume in terms of the variable z/, 

V 



(2.6) 



where 7 is the Barbero-Immirzi parameter, Ip = \jGh is the Planck length (we take c = 1), and e — ±1 determines 
the orientation of the physical triad relative to the fiducial (co-)triad determining qab (= '^a'^l^ij) ~ see Refs. 
[IIEIEIE]. Thus —00 < V < +00. Note that v is dimensionful. For comparison to other work, note that v = \ - 
where v is the dimensionless volume variable of Refs. [5j[6], and 

\^yfK-lp (2.7a) 
4^3777 • Ip. (2.7b) 
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Here A • is the "area gap" of loop quantum gravity! 

Remarkably, when the physical model given by Eq. (2?2) is loop-quantized in these variables, the classical "harmonic" 
gauge choice N{t) — a{t)^ leads to an exactly solvable quantum theory]^ referred to as "sLQC" (for "solvable LQC") 
[Uin]. One finds that physical states 4* (1^,0) may be chosen to be "positive frequency" solutions to the quantum 
constraint, 



(2., 



where the positive, self-adjoint "evolution operator" (the quantized gravitational constraint) O is given in the h'- 
representation by a second-order difference operatorj^ 



(e*)(z/, 0) = {v/|j^(i^ + 4A)||t/ + 2A|^'(iy + 4A, (j)) - 2iyH{i^, dp) + ^\v{v ~ '^X)\\v - 2X\^{v - 4A, </.)} . (2.9) 

Solutions to the full quantum constraint (C = — [9^ -f 6]) therefore decompose into disjoint sectors with support on 
the e-lattices given by = 4A?t. -I- e, where e £ [0,4A) iSlE]- In order not to exclude the classical singularity at = 
from the start, we work exclusively on the lattice e = 0, so that in this quantum cosmological model, the volume is 

discrete: 



V — 4An, n Cz "Zi. 
Group averaging yields the physical inner product 

{^\^)= ^ vl/(j.,0„)*<i>(j.,^„) 



(2.10) 



(2.11) 



for some fiducial (but irrelevant) According to Eq. (2.8), states at different values of the scalar field (p may be 
mapped onto one another by the unitary evolution 



(2.12) 



It is natural therefore - though not essential [S] - to regard the scalar field as an emergent physical "clock" 
or "internal time" in which states evolve in this model. Eq. (2.121 shows that the inner product of Eq. (2.11) is 
independent of the choice of (jjo, and is therefore preserved under evolution from one (/(-"slice" to another. 

Finally, we note that in the absence of fermions, the action, dynamics, and other physics of the model are insensitive 
to the orientation of the physical triads [21 [51 IHl [TS] . We may therefore restrict attention to the volume-symmetric 
sector of the theory in which 



*(j^,0) *(-z/,, 



(2.13) 



Many further details concerning the quantization of this model and its observables may be found in Refs. [D [21 [3 



MM- 



C. Observables 



The basic variables in this representation are the scalar field and the volume v. Employing (j) as an internal time, 
the primary operators of interest are the volume, which acts as a multiplication operator, 

= (2.14a) 
V^{i,,^) = 27rjll\i^\^{i^,<j>), (2.14b) 



^ Note that in earlier work in loop quantum cosmology A was given as 2\/3n^. However, in Ref. [J it was sliown tliat A sliould be 
taken to have twice that value in homogeneous models. Since the volume eigenvalues were given in terms of the area gap this difference 
does not intrude unduly into those prior results. The relation u = X ■ v holds so long as one employs the same area gap consistently 
throughout. 

^ The choice of the harmonic gauge leads to the exact solvability of the quantum theory because it eliminates inverse factors of a in the 

Hamiltonian constraint; cf. Eq. | |2.3[ l. 
^ This expression is different from what is found in Refs. [S] [B] because we are using states that carry an additional factor of y'A/|j^| 

relative to those states in order to simplify the form of the inner product, Eq. ( |2.11| . Compare, for example, Refs. [TJE]- 
® See footnote [5] 
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and the scalar momentum p^, 

P0 ^{ly, (f)) = -ihd^^{v, (f)) (2.15a) 
= ;i\/e^^'(j/». (2.15b) 

(In this paper we will not have need of the (exponential of the) momentmn b conjugate to [51[S].) As in the classical 
theory, the scalar momentum is a constant of the motion - it obviously commutes with the effective "dynamics" 
given by ^/Q ~ and is therefore a Dirac observable. The volume v is not, but the corresponding "relational" observable 
giving the volume at a fixed value (p* of the internal time <j) is. Defining 

[/(0) = e*^"^, (2.16) 
the "Hcisenberg" operator i>|0. (0) acting on states at (f) is given by 

^)U.(0) = [/(0*-0)^z>C/(0*-<^), (2.17) 
so that, for example, the physical volume V — ^ttjI^ of the cell V at </)* is given by the operator 

V\^, {P) ^{v, 0) = 2n-ille'^"^^-^'^\v\^{v, 0*). (2.18) 
It is straightforward to verify that p^ and V\^-{(j3) commute with U{4>), and are therefore Dirac observables. 

III. EIGENFUNCTIONS OF THE EVOLUTION OPERATOR 

General physical states ^'(j^, (j)) may be readily expressed in terms of the cigenfunctions of the dynamical evolution 
operator Q - the gravitational part of the Hamiltonian constraint - given by 

Q^ek{iy) = ujlekiv), (3.1) 

where 



uJk^Vl27TG\k\ (3.2a) 
= K\k\, (3.2b) 

and — oo < fc < oo is a dimensionless number labelling the 2-fold degenerate eigenvalues. Restricting to the symmetric 



lattice v — 4An and physical states which satisfy Eq. (2.13), we usually choose to work with a symmetric basis of 



cigenfunctions e^if\i') which satisfy e^^\iy) — e\^\—iy). In terms of these physical states may be expressed simply as 

/ + 00 
dk^{k)ei"\iy)e"^'''''. (3.3) 
-OO 

For normalized states, J2iy^AXn — 

f + 00 



d/c|*(fc)|^ 1. (3.4) 

Explicit analytic expressions for the cigenfunctions of Q have recently been found [7j |8] . The symmetric cigenfunc- 
tions will eventually be expressed in terms of the primitive cigenfunctions [7] 

eo{'^)^^Q.M (3.5a) 

efc(i^) = A(fc)-^/5d^ ^ d6e-*'^e*'=i"(*""^) (fc ^ 0), (3.5b) 
where A(k) is a normalization factor which for consistency will always be chosen to be 

A{k) ^ ^ (3.6) 
^y Airk sinh(7rfc) 
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The functions efe(i^) have support on both positive and negative v and are not symmetric in v. It is convenient to 
seek hnear combinations e^(y) of e^iy) and e-k(y) which have support only for ^ 0. The correct combinations 
turn out to be [7] 



Clearly 'Y^y^^iy')* e^,(y) — 0. The choice of A(fc) in Eq. (3.6| corresponds to the normalization 

E e±M*e±(z.)=5(^)(fc,A:'), 



i/=4An 



where J*^"^ (fc, fc') is the symmetric delta distribution 



(3.7) 



(3.8) 



S^^){k, k') = - {6{k, k') + 5{k, -fc')} . 



(3.9) 



(In contrast to Ref. [7], we choose to work with the full range offc, — oo < k < oo. This leads to the second 
delta function appearing in Eq. (3.8) relative to Eq. (C13) of that reference. Since as we will see these functions 
are symmetric in k, the two approaches are of course equivalent, but do lead to some differences in choices of 
normalization.) 



The e J (v) can be given explicitly as [7] 



7r|iy| 
A 



e^(i/) = A{k)\j ^I{k, ±!./4A) 



A{k)J'^I{k,±n), 



recalling ly = 4An. Here /(fc, n) = for n < 0, and for n > is given bjj^ 



2n 



1 



2n-l 



I(k, n) — iky - , s , 

^ ' ' ^ to! 2n-m)! 

m=0 ^ ' 1 = 1 



J]^ (ik + m-l) 



= —ik 



r{2n-ik) 



r(l + 2n)r(l-ifc) 



—T 2Fi{ik, -2n; 1 - 2n + ik; -1), 



(3.10a) 
(3.10b) 

(3.11a) 
(3.11b) 



where the second form follows from the first by simple manipulations of the definition of the hypergeometric function 
2Fi(a,6;c;z) [16]. 



We will discuss the properties of I{k,n) in detail later. For now, note from Eq. (3.10) that 



and from Eq. (3.7) that 



e±,(^)=e±(z.). 



Given the symmetry relation Eq. (3.12) it is clear that the symmetric eigenfunctions e^f^v) are finall 

1 



V2 



{e+(^) + en^)} 



A{k)\l^ I{k,H/4X) 



2|A:sinh(7rA;)| 



I{k,\n\). 



(3.12) 
(3.13) 

(3.14a) 
(3.14b) 
(3.14c) 



Compare Ref. lZj, Eq. (C8) and Ref. g], Eq. (B4). 
* Notice that restricting to the volume-symmetric sector of the theory, Eq. ( |2.13| , on the e = lattice lifts the 2-fold degeneracy of the 
eigenvalues oj^ to a single eigenvector for each k [6]. 
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The following symmetry properties may be verifiedj^ 



The e^'^(i/) with A{k) chose n as in Eq. i3.6\ then satisfy the completeness relations 



+00 

— 00 



(3.15a) 
(3.15b) 
(3.15c) 

(3.16a) 
(3.16b) 



where the symmetric Kronecker delta is defined analogously to Eq. (3.9 1. (The domains of these expressions are 
understood to be even functions of k and u, respectively.) The symmetrized deltas arise because the functions e^^\v) 
are symmetric in both v and k. 



An expression for e^,'''' {v) we will find useful later is 



cosh(7rA:/2) 
71 



{ek{v) + e_fe(j/)} . 



(3.17) 



which follows from Eqs. (3.14a) and (3.7). As a useful aside, note it is easy to see from Eq. (3.5) that e'^\h')* — 
ei}^{—v). Additionally, the change of variable b' — —b + tt/A in Eq. (3.5) - remembering v — AXn - reveals that 

eki-iy)=e_kii^). (3.18) 



Thus efc(j^)* = e-ki—v) — efc(i^), and both efe(i^) and e"^\i') are therefore real. 

This completes the catalog of properties of the eigenfunctions we will require. We now describe the behavior of 
the functions e^^\v) we seek to explain in the sequel. The results of the analysis will confirm and complement the 
understanding of earlier numerical and analytical work arrived at prior to the discovery of the exact solutions for this 
model. 




(a) n = 20 




(b) n = 200 



FIG. 2. Plotofe<^-'>(i/ = 4An) as a function of fc > for n = 20 and n = 200. Regarded as a function of fc, e)^ {v) is symmetric in 
k and always exhibits exactly n — 1 nodes in addition to the node at fc = 0. The largest zero and last maximum of e^^\v) always 
appears at a value < 2|n|, after which e]^\u) is exponentially damped in k. This ultraviolet cutoff at |fc| = 2\n\ = |i'/2A| 

in the eigenfunctions will be explained analytically in the sequel. 

Plots of e^^\v) are shown in Figs. 2 3 and 4 Two behaviors are clearly evident in these plots. First, the dynamical 



eigenfunctions are exponentially damped for^fcl > 2\n\ = |j^/2A| (Figs. 2|3). This is the ultraviolet momentum space 



^ Eq. | |3.15a^ follows from Eq. l |3.14| l. Eq. | |3.15b| l follows from Eq. l |3.13| l. Eq. H .15c| follows from a study of I{k, \n\) or the argument 
to follow. 
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cutoff in the eigenfunctions described in the introduction. The cutoff can be understood as a consequence of the 
fundamental discreteness of volume in these quantum theories. Second, for \n\ > \k\, the eigenfunctions e^i^\i^) settle 
quickly into a decaying sinusoidal oscillation (Fig. |4]). We will see that this oscillation corresponds to a specific 
symmetric superposition of the eigenfunctions for the Wheeler-DeWitt quantization of the same physical model. As 
a consequence, generic quantum states in this loop quantum cosmology will evolve to a symmetric superposition of 
an expanding and a collapsing Wheeler-DeWitt universe. 



0.3-— *♦ 



100 




FIG. 3. Plot of the functions el!'\iy = 4An) in the (fc,n) plane for < n < 100 and |fc| < 100. They are symmetric in both k 
and n. The volume variable f = 4An is fundamentally discrete; only discrete values of the continuous variable k are plotted 
for simple reasons of visual clarity. The exponential ultraviolet cutoff along the lines |fc| = 2\n\ = is clearly evident. 

The dynamical eigenfunctions e^i^\i/) may therefore be regarded to an excellent approximation as having support only in the 
"wedge" |fc| < 2\n\. It is this feature of the eigenfunctions that is ultimately responsible for the existence of a universal upper 
bound to the matter density. 



A. Asymptotics 



It is evident from Fig.jsjthat, to an excellent approximation, the dynamical eigenfunctions el^\iy) may be regarded 
as having support only in the wedge |fc| < 2\n\ in the {k,n) plane. We now seek to explain this behavior based on an 
analysis of th e exac t sol utions , Eq. ( |3.14 1, as well as the large volume limit and other features visible in Figs. 2]|4 



From Eqs. (3.14) and (3.111, el^\iy) is given by ^/\n\/2\k ■ sinh(7rA:)| times the polynomial /(fc, \n\). Indeed, exam- 



ination of /(fc, regarded as a polynomial in k, shows that it is an even polynomial in k with no constan t term . 



whose terms alternate in sign. Thus we see nnmediately that 6^*^(0) = and e^^^(j/) = for 9^ (cf. Eq. (|3.5a|)). 



Examination of plots of I{k, |n|) (as in Fig. [2]) shows that it always exhibits the maximum number of roots possible 
(2n — 1) for such a polynomial; the alternating signs of the coefficients lead to the oscillations. 
Writing the product in Eq. (3.11a) as a "falling factorial" ^16j, it is possible to arrive at an explicit expression for 
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(a) fc = 10 (b) A: = 100 



FIG. 4. Plot of e^^\v = 4An) as a function of n > for fc = 10 and k = 100. Regarded as a function of n G Z, the ej^ {u) are 
symmetric in n and have support only on the lattice v — 4An; the points are connected for visual clarity only. Note e^°^(0) = 
for fc 7^ 0. The functions e'^i^\i' ~ 4An) for fixed k decay rapidly to essentially zero for |n| < |fc|/2, the ultraviolet cutoff in the 
eigenfunctions also visible in Figs. |4]|3] (The cutoff is not as sharp viewed on slices of constant k as it is on slices of constant 
n, on which the cutoff is truly exponential.) For |n| > |fc|/2, they settle rapidly into a regular decaying oscillation. This latter 
behavior corresponds precisely to a symmetric superposition of Wheeler-DeWitt eigenfunctions to be elaborated in the sequel. 



/(fc, |n|) which is useful for some computations 



\n\ 

J(fc,|n|)=^a(|n|,j)fc'^ (3.19) 



where the coefficients a(|rt|, j) are given by 



.(„.;) . (-1)^ E ,(2„ - M - 1) . f ^ J . E i!;!^ ) . (3.20) 



l=2j 



Here the s{p,q) denote the (signed) Stirhng numbers of the first kindj^ It should be noted that the factor to the 
right of (—1)-' is always positive, leading to the alternating signs of these coefficients. 
For large I{k, \n\) is dominated by fc^l"l, and therefore 

e\!^ {v = 4An) - ^ • fc^l"! (3.21a) 

i/fc sinh(7rfc) 

- •6"'^'=/^ (3.21b) 

and the decay of the eigenfunctions is indeed exponential in fc past the largest root of e'^^iv). 



Steepest descents 



To identify the value of k at which the decay of the symmetric eigenfunctions e^^^(i^) sets in requires a bit more 



work. Recall from Eq. (3.17) that the el^\h') may be expressed in terms of the primitive eigenfunctions ek(v) given 



by Eq. (3.5 1. The integral in this equation is of the form 



I(fc, v) 



db e^^^^'*-''") 



(3.22) 



One must be careful, however. The alternating signs of the coefRcients and the large powers of k appearing in the expression for I{k, |n|) 
at large volume (n) can quickly lead to numerical instabilities. 

The rising and falling factorials are defined by x" = r[r=(/(^ ~^ *) ^'^'^ ^" ~ Yl.7^o(^ ~ '^'^^ signed Stirling numbers of the first kind 
are then defined by x— = X]ILo s{n,i)x^. 
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where 



fib, k,v) = k- ln(tan y ) - y . 



(3.23) 



This was the form from which the exact expression Eq. (3.101 was extracted in Ref. [7]. We however wish to evaluate 
this integral in the limits of large |fc| and \v\. While Eq. (3.22) is not quite of the same form for which the steepest 



descents approximation is normally discussed - a single large parameter multiplying an overall phase - the same 
arguments for the validity of the approximation apply. In regions where k, v) is large, the integrand oscillates 
rapidly and contributions from neighboring values of h cancel one another. The dominant contributions to i^), 
therefore, come from regions close to the stationary points of fc, v) where / changes only slowly with h and the 
cancellations are not strong. This is the usual steepest-descents approximation, and in general one has !17] 



2tt 



\n^^)\ 



i^signif" {zi)) 



(3.24) 



where the Zi locate the stationary points f'{zi) = along the relevant contour. 

It is clear that when |fc| or are large, f{b,k,i>) can become large, suppressing the value oi T{k,iy), and so we 
seek the stationary points of /. 

First observe that /(&, k, v) diverges at 6 = and b = tt/A, so there is no contribution to T{k, v) from the endpoints 
of the integration due to the rapid oscillation of the integrand there. Next, one finds 



dh 



{b,k,u) 



Xk 
sin Xb 



(3.25) 



and 



5V 

a&2 



{b,k,v) 



X^kcosXb 
sin^ Xb 



(3.26) 



The stationary points therefore satisfy 



sin A6 = 



2Afc 



(3.27) 



When solutions exist there are two roots bi and &2, given in the limit 3> A|fc| by 



bi 
b2 



2k 

A " 



2k 

V 



(3.28a) 
(3.28b) 



In this limit 



f{bi,k,i^) w -fc 
f{b2,k,:y)^+k 



V 

Xk 

V 

Xk 



TTl/ 

2A' 



(3.29a) 
(3.29b) 



and 



— ib„k,.)^--, 



(3.30a) 
(3.30b) 



We now piece together these results to study the asymptotic limits of X(fc, v) and consequently e^^"^ (v) 
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2.. Wheeler-DeWitt limit 



We will begin by considering the case in which k and v are both positive, and return to the other possibilities 
shortly. When > Afc > 0, from Eq. p. 241) we find 



I{k,v) 



27r 



if{bi) ^ifsign{f"{bi)) ^ 



27r 



2^/^cos(fc 



l/"(&2)| 





+ 4) 







(3.31a) 

(3.31b) 
(3.31c) 



where to get to the last line we recall v = AXn, so the final exponential factor in Eq. (3.31b) is unity. 
In the case where k and v are both negative, the same results obtain, but now 



f{h,k,:^)^+\k\ 
f{b2,k,,.)^-\k\ 



2A 



sign{J"{h,)) = +, 



again leading to a cosine, but with k — > |fc|. Thus, from Eqs. (3.5), (3.6), and (3.31c), we fincp] 

■0' 



2A 



yj\ sinh(7rfc)| 



cos I I A: I In 



a(|fc|) ) when 



\v\ > A|fc| 
v-k>Q 



where 



t{k) = fc(l - Infc) + 



(3.32a) 
(3.32b) 



(3.33) 



(3.34) 



We have yet to consider the case where k and v are opposite in sign. In this case note that since < 6 < tt/A, there 
are no solutions to Eq. (3.27), and /(6, fc, v) has no stationary points in the domain of integration. Thus I(fc, j/), and 
hence efc(^), are strongly suppressed by the rapid oscillations of the integrand when k and v are opposite in sign. 



We note from Eq. (3.17) that the symmetric eigcnfunctions e'^\v) are a linear combination of ek{v) and e_fc(i^). The 



functions e-fc (t^) a re, mutatis mutandis as above, strongly suppressed when v and k have the same sign, and assume 
the limit Eq. ( 3.33 ) when k and v are opposite in sign. The \v\ 3> A|fc| limit of e^,''"' {v) will therefore pick up precisely one 

no matter the signs of k and v. Observing that cosh(7rfc/2)/Aysinh(7r|fc|) w 1/^/2 
, we arrive finally at 



contribution of the form of Eq. ( 3.33 
for even very modest values of k > 



As) 



(^) 



/ 2A 
tt\v\ 



(|fc| 



In 



-a(|fc|) 



\v\ > A|fc|. 



(3.35) 



Fig. [5] shows that the exact eigenf unctions settle down to this asymptotic form very quickly. We will employ Eq. 
( 3.35[ ) to study in Sec. IV the large volume limit of flat scalar loop quantum universes. 

The asymptotic expression Eq. (3.35 ) was, in effect, arrived at on the basis of analytical and numerical cons iderat ions 
in Ref. with a numerically motivated fit for the phase a(A;)j^ In Ref. ,11] an expression equivalent to Eq. (3.35 ) was 
derived from a careful analysis of the asymptotic limit of solutions to the constraint equation, including an expression 
for the phase a{k) equivalent to Eq. (3.34). (See Ref. [15] for a related analysis of this limit.) Here we have instead 
derived this asymptotic form from the exact eigcnfunctions, explicitly confirming these prior analyses with the exact 
solutions for the model. 



We could also have arrived at the absolute value signs simply by observing from Eq. (3.181 that e_k{—v) = ef,{v). 
Consult footnote [5] concerning the factor leading to the decay of the eigenfunctions with increasing volume. 
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FIG. 5. Plot as a function of n of both the dynamical eigenstate e^^\v 



The volume variable v = 4An is fundamentally discrete; the values of e^"' (n ) 
by a dashed blue line for visual clarity. The solid red curve is the corresponding asymptotic form. The rapid convergence to 
the asymptotic form for |n| ^ on the lattice f — 4Xn is clear. Note this asymptotic form c orres ponds to the particular 
superposition of eigenstates of the Wheeler-DeWitt quantization of the same model given by Eq. (4.5 1. The rapid oscillations 
visible at small volume are the correct physical behavior of the Wheeler-DeWitt states, and are ultimately responsible for the 
fact that these models are singular. See Refs. |14lll8j for further discussion. 



4An) and the asymptotic form Eq. (3.351 for k = 30. 
are marked with blue x's; the points are connected 



3. Ultraviolet Cutoff 



Figures bis clearly exhibit the exponential ultraviolet cutoff in the eigenfunctions e^^\u) for values of |fc| > 2\n\ 



v/2\\. We know already from Eq. (3.21) that an exponential decay will eventually set in. The only question is 



at 



what value of k does that occur? We have, in fact, already seen the origin of this cutoff and its value. Eq. (3.27) 



shows that f(b, fc, \v\) has no stationary points when |2Afc/z^| > 1. In other words, I(A:, v), hence efc(i^) and ej^, ( 
are strongly suppressed unless 



|fc| 



< 



J/ 

2A 

= 2|n| 



(3.36a) 
(3.36b) 



This cutoff - in particular, its linear scaling with volume - may be understood physically as a consequence of the 
underlying discreteness of the quantum geometry. States with wave numbers |fc| > 2|7z| (i.e. wavelengths shorter than 
the scale set by |A/i^|) are not supported. Alternately, it may be viewed as the manifestation in the eigenfunctions of 
the "quantum repulsion" generated by quantum geometry at volumes smaller than the wave number. 



4- Small volume limit 

The same argument shows that the eigenstates will, equivalently, decay rapidly for small volume, when |n| ^ |fc|/2, 

tells us the decay in e^,"^ (v) as a function of k is exponential. The precise functional 
form of the decay as a function of n is less evident, but the figures show it is also quick. 

At this point a comment may be in order. It is tempting to study this question by regarding e^i^\i') as a function 
of a continuous variable v. However, plotting the exact expressions for el^^v) for continuous values of v on top of 

the values for v = AXn should quickly disabuse one of the notion that there is a simple sense in which e^.''^ (v) is well 
approximated by its naive continuation to the continuum. In fact, as discussed in detail in Ref. [5], the convergence 
to the Wheeler-DeWitt theory in the continuum is not uniform, and must be extracted with some care in the limit 
the "area gap" set by A - fixed in loop quantum cosmology to the value of Eq. ( |2.7[ ) - tends to 0. 

As noted in the Introduction, and as is clearly evident in Fig. [3] the ultraviolet cutoff in momentum space is the 
"same" cutoff as the rapid decay at small volume as a function of volume that has long been known in loop quantum 
cosmology based on numerical solutions for the eigenfunctions [51 H] . It was also known numerically that the onset of 



as is clear in Fig. wj Eq. (3.21 
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this decay was proportional to the eigenvalue Wfc. (See e.g. [Ulin]-) What is new in the present work, facilitated by 
the change in perspective to consideration of the behavior of the eigenfunctions as functions of the continuous variable 
k, is an analytic understanding of the linear cutoff grounded in a study of the model's exact solutions, and its specific 
relation to the critical density. 



B. Representation of operators 



As noted above in Eq. (2.13), we have restricted attention to the volume-symmetric sector of the theory. This is 
only possible because the physical operators preserve the symmetry of the quantum states. 
From Eq. (2.9 1, the matrix elements of Q in the volume basis may be expressed as 

{iy\QW') = 12TrG^\n-n'\\n + n'\ ■ |,5„,„, - ^ [Sn,n'+i + Sn.n'-i]j , (3.37) 

where v = 4An and = 4An'. Note these matrix elements satisfy the following properties: 

{iy\e\iy) = (-i/|e|-i^) (3.38a) 
= W\&WT (3.38b) 
= (i^'|e|i/). (3.38c) 

Owing to these relations, the operator 9 preserves the subspaces "Hp^y^ and 'H'^^y^ of states that are even and odd in 
1^, so that p('^)ep('^) = 0, where P'^"'^ and P(") are the corresponding projections. (In other words, 6 commutes with 
the parity operator 11^ = P^''' — P'"' 011].) On the symmetric subspace "Hp^^g to which we have restricted ourselves, 
6|^(s) = p(s)Qp(«) = qI") (and correspondingly p^*'' = hV 0(*)) may be decomposed in terms of the symmetric 

phys 

basis of eigenstates 

Q^'^ J dke\k'^'^){k'^'^\, (3.39) 
where e^i^\iy) = (i^|fc(^)). The matrix elements of O'^''^ in the volume representation are related to those of 8 by 

(HQ^^V) - I {(Hei^') + (Hei-^')} (3.40a) 

= {iy\Qi^^\-v'). (3.40b) 

The actions of O and O^'*^ on Hpjjyg are of course completely equivalent. Since p"^ — h?Q on "Hphys, these expressions 
give the matrix elements of on "Hp^y^ as well. 

The matrix elements of = h^/Q are more complex. These are given in terms of derivatives of a generating 
function in Appendix C of Ref. |7 . Explicit expressions for the physical observables in another representation are 

also given in Ref. [3]. Here we note that on "Hp^yg we may employ the e^^\v) to calculate ff^'^ explicitly in the volume 
representation. Indeed, the polynomial solution Eq. (3.19) for I{k,n) makes it a straightforward matter to evaluate 
these matrix elements. The result is (with v — 4An and u' — 4Am) 

r dk \k\ ei'\,y)ei'\,.r = VWM E E «(l"l'^>(l"^M)x5ji^^^^(^r(2(j +l) + l)C(2(j + /) + 1) 

j=i 1=1 

(3.41a) 

l"l I™! 2 
« y/\n-m\ E E ai\nlj)ai\mll)-^^j^ri2{j + I) + 1), (3.41b) 

j=i 1=1 

where C(z) is the Riemann zeta- function. This expression gives the {ly, v') matrix elements of p^^^ /Hk, or equivalently 
a/oW/k. 

We observe from Eq. (3.37) that is nearly diagonal in the volume representation, with only the n' — n, n zfc 1 
elements not exactly zero. The same is therefore true of p\. Direct numerical evaluation of the expression Eq. (3.41 ) 
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T"^ — 




FIG. 6. Plot of the norm of the scalar momentum matrix element overlap integral appearing in Eq. (3.451 normalized by 



the estimated maximum value }^\v/2\\ = \n\ of this integral on the diagonal \v\ — taking = 4Xn\ < \v' = 4Am 



\ J^^dk\k\e^^\v)e]^\v'Y\/\n\. The integrals have been calculated numerically from the exact eigenfunctions over the range 



< \n\ < 50 and < m < 50. According to the upper bound expressed in Eq. (3.441, this normalized matrix element is 
bounded above by one, and as argued is strongly suppressed off the diagonal \v\ = \v'\. In effect, these plots show that 
approximately commutes with since it is nearly diagonal in the i^-representation. This is essentially the reason the "moral" 
argument expressed in Eq. (5.81 for the existence of a critical density in this model yields the correct result. 



reveals that - and therefore \/Q - are also nearly diagonal in the volume representation, with only the n' = n, n± 1 
matrix elements significantly different from zero. (See Fig. [6)) In this case, however, the ofF-diagonal elements of 
are merely very small, rather than precisely zero. 

The values of these matrix elements can be understood as a consequence of the ultraviolet cutoff, Eq. (3.361. Indeed, 
the exponential cutoff < |i^/2A| implies that the diagonal matrix elements are bounded 



dk\k\e'^^\v)e^^\±vr 



< 



V 

2A 



(3.42) 



(The 1/2 is a consequence of the symmetric normalization of the eigenfunctions, Eq. (3.16 1.) This bound on / k 

may be compared with that set by the exact expression for O, Eq. (3.371. From Eq. (3.40), 

{iy\e^'^\u) = l{v\e\u) (3.43a) 



K 



2A 



(3.43b) 



Numerical integration of the exact eigenfunctions shows that this rough estimate of the upper bound typically overestimates the actual 
value of the integral by around 60% due to the (suppressed) contribution to the normalization integral from the increased amplitude 



near |A:| ~ W /2\\. See Fig.je] 
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As it is always the case that (A^) > {Ay, we see that a strict bound on the diagonal matrix elements of VQ^/k is 



dk\k\e'^^\j^)el'>i±iy 



1 

< 

- 2 



2A 



(3.44) 



in agreement with the bound inferred from the UV cutoff. 

The off-diagonal elements may be bounded in a similar manner. For simplicity assume |^'| < {i^'l. The exponential 
UV cutoff effectively restricts the range of integration to |A;| < |i^/2A|. The Cauchy-Schwarz inequality then gives 



dk\k\e^^\iy)ei^\:.' 



< 



V 

2A 




dk\et\y)?- 



~\vl2\\ 



dk\e^:\v'W. 



(3.45) 



Again, because the e"^\i') are symmetrically normalized, the value of the first square root is essentially l/\/2. As 

for the second, we note from Fig. [2] that the e^*'' [v) execute approximately uniform amplitude oscillations, growing 
slowly with increasing k with a short lived increase before the exponential cutoff sets in at |fc| = |i^/2A|. Therefore, 
for \v\ < 1 1/' I we may estimate that at most 



k/2A| 



k/2A| 



(3.46) 



showing that the cutoff alone implies that the off-diagonal terms are suppressed relative to the diagonal terms. 
Interference effects only reduce their values further; Fig. [6] shows that except for the n = n' ±1 elements - as with Q 
itself - this suppression is dramatic. 

As a shorthand to express these bounds, we can say that and - and therefore |£L - approximately commute, 
in the sense that {v\p^\i/') is approximately diagonal. (See Fig. |6]) We will see in Sec. |v]that this helps explain why 
the matter density in these models remains bounded even though the spectrum of the scalar momentum p^ is not 
itself bounded. 



IV. LARGE VOLUME LIMIT OF LOOP QUANTUM STATES 



In Eq. (3.351 we have exhibited the large volume (more precisely, ^ A|fc|) limit of the basis e^, {v) of symmetric 
states of nat scalar loop quantum cosmology. We extracted this limit from the exact solution for the model's eigen- 
functions, essentially confirming prior numerical and analytical work. In this section we relate these states to the 
eigenstates in a Wheeler-DeWitt quantization of the same physical model. 

A complete, rigorous Hilbert space quantization of a flat Friedmann-Lemaitre- Robertson- Walker cosmology sourced 
by a massless minimally coupled scalar field has been given in Refs. [U [6] and compared to its loop quantization in 
detail in Ref . [3] . It is known rigorously that states in the Wheeler-DeWitt quantization are generically singular just 
as they are in the classical theory in the sense that all states assume arbitrarily small volume (equivalently, large 
density) at some point in their cosmic evolution in "internal time" (j) HI] . 

The classical solutions are given in Eq. (2.4), corresponding to disjoint expanding and contracting branches which 
either begin or end in the classical singularity at = 0. Solutions to the Wheeler-DeWitt quantum theory similarly 
divide into disjoint expanding and contracting branches, and as noted, are singular in the same way. 

The Wheeler-DeWitt version of the quantum constraint i^^ [6] 



127rG 



WdW WdW 



(4.1a) 
(4.1b) 



Attention may again be restricted to symmetric (Eq. (2.13)), positive frequency solutions in the sense of Eq. ( |2.8 
The symmetric eigenstates of satisfying Eqs. (3.1)-(3.2) are 



1 



ik In 



(4.2) 



See footnote [5] 
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and are orthonormal (distributionally normalized to S{k, k')) in the inner product 



WdW I ^WdW 



resulting from group averaging. Physical states may then be expressed as 

/+00 
-oo 



(4.3) 



(4.4a) 



1 

-j^ /"OO 



(4.4b) 
(4.4c) 



^4^70 
EE ^'™"(z/,0) + ^'™(zy. 

The orthogonal sectors of "right-moving" (in a plot of ((> vs. v) and "left-moving" states clearly correspond to the 
expanding and contracting branches of the classical solutions, Eq. (2.4|. A priori, note that ^^"^^(fc) = '^^'^^{k) {k < 
0) and *^''*(fc) = ^^-"^{k) (fc > 0) need not be in any way related in the Wheeler-DeWitt theory. 

We now show that generic states in the loop quantized theory decompose into symmetric superpositions of expanding 
and collapsing (right- and left-moving) Wheeler-DeWitt universes at large volume. This follows simply from Eq. (3.35 1, 
which may be written 



(4.5) 



where z = y/2 for the normalization of Eq. (4.2 1 appropriate to the range —oo < v < oo of Eq. (^4.3M^^ and the -\/A 



is present because the e'j^\u) are dimensionless, whereas the e^''^(i/) are not. The factor of zy/\ can be understood 
as arising from the difference between normalization of the Wheeler-DeWitt eigenfunctions on the continuous range 
—oo < V < oo vs. the normalization of loop quantum states on an infinite lattice with step-size 4A. (Compare 
Appendix B of Ref. [H].) 



The relationship expressed in Eq. (4.5) has long been known in loop quantum cosmology on the basis of both 
analytic and numerical arguments. See, for example, Eq. (5.3) of Ref. [6], as well as Eq. (3.1) of Ref. [11] - which also 
contains a careful analysis of the convergence properties of this limit - among many others. Here we have confirmed 
that the asymptotic behavior of th e exact solutions agrees precisely with these earlier arguments. 



We know that the limit Eq. (4.5 1 is valid when \v\ ^ 2A|A;|, and more generally, that e'^\v) has support only in the 



wedge |fc| < \v\/2\. E q. (|4.5[ ) therefore holds inside this wedge of support but clearly breaks down near its boundary 
|fc| = 2|n|. From Eq. (3.3), quite generally 



dk ^{k)e^^\v)e"^^^ 

-OO 

kl/2A 



(4.6a) 



dk ^{k)e^^\v)e''^'''l'. (4.6b) 

'-kl/2A 

It is noted in Ref. [TTJ that one must take care to draw conclusions concerning the asymptotic behavior of states in 
sLQC based on that of the eigenfunctions because the convergence of the sLQC basis e'^\v) to that of the Wheeler- 
DeWitt theory is not uniform in k. Nonetheless, we argue that for a wide class of quantum states, there will be a 
well-defined region depending on the state in which this approximation will hold for that state. 



Specifically, replacement of e'^\v) in the expression Eq. (4.6) with its asymptotic form Eq. (4.5) will be valid for 
values of the volume (significantly larger than that) for which the Fourier transform ^'(fc) does not have significant 
support outside the wedge at that volume. Quantum states are normalized, so we know that '^(k) is square- integrable. 
Because functions of compact support are dense in L^(M), there is a dense set of states for which, for every state 
^(z^, (/)) in this set, there is some value of |fc| - call it k^^ - outside of which ^(fc) has no support. Therefore, for a 
dense set of states in the quantum theory the replacement Eq. (4.5) will be a good approximation for \v\ 3> 2A|fc^|. 



^® Since states 0) are symmetric in v, one is free to restrict instead to positive u only, in which case the v^Stt in Eq. i 4.2 i should 

be V27r and z = 1. 
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It is worth emphasizing, therefore, that the domain of applicabihty of the large-volume approximation is dependent 
upon the quantum state through the support of "^(k). 

For states satisfying this condition and within that domain of applicability, we may write 



\v\l2\ 



\u\/2\ 



(4.7) 



It will be seen short ly th a t the first t erm co rresponds to a contracting universe, and the second, expanding. To begin, 
we n ote from Eqs. (3.2|, (3.3), and (3.15b) that 4'(fc) is even, \t'(— fc) = ^'(fc). Consider the first term alone. As Eq. 



(4.7) applies only for values of the volume for which '^{k) has negligible support for |fc| > |j/|/2A, we may extend the 
range of fc-integration to — cxd < k < oo. By separating the integral dk ~ dk + dk and making the change 
of variable k' — —k in the first, one quickly finds 



(4.8a) 
(4.8b) 



As in the Wheeler-DeWitt case, Eq. (4.4), ^'^(i^, 0) clearly corresponds to a contracting quantmn universe, with 

(4.9) 



equivalent Wheeler-DeWitt Fourier transform 



^WdW 



(k) = 5'(fc)e' 



4\k\) 



fc > 0. 



In an exactly similar way, the second term in Eq. (4.7) is 



^R{iy, (j)) = 2zVX / dk *(/fc)e~*" 

J —OO 



WdW 



''.K\k\cl) 



which describes an expanding universe with equivalent Wheeler-DeWitt Fourier transform 

^-^"-"(k) = §(fc)e-*"(l'^'l) k < 0. 



(4.10) 



(4.11) 



Therefore, for a dense set of quantum states and within the domain of applicability of the large volume approximation 
1 1^1 ^ 2A|fc*| for the state ^'(z/, (j)), we may always write 



vl/(z.. 



*i?(i^, 



*l(j^,< 



large 
volume 



(4.12) 



Unlike the Wheeler-DeWitt case, however, 5'^ and '^ji are not independent. In fact, owing to the symmetry of '^{k) 
in the loop quantum case, the equivalent Wheeler-DeWitt Fourier transforms are essentially the same, having equal 
modulus |^l(— fc)| = |^fl(fc)| and a fixed phase relation given by exp(iQ;(|A:|)) between them. (Note that expressions 
equivalent to Eqs. (4.9) and (4.11) may also be found in Ref. [IT]. These relations are central to the "scattering" 
picture of loop quantum cosmology developed in that reference.) 



Thus, we have shown from the exact solution that in the sense given by Eq. (4.12), at sufficiently large volume a 



dense set of states in flat scalar loop quantum cosmology may be written as symmetric superpositions of expanding 
and contracting universes. This is an essential feature of loop quantum cosmology, deeply connected with the fact 
that these cosmologies are non-singular [T5] - all states "bounce" with a finite maximum matter density. 

It was observed long ago from numerical solutions that semiclassical states "bounce" symmetrically, including the 
dispersions of these states jSl |6] . Analytic bounds on the dispersions of all states in this model have been proved in 
Refs. [m [12] , in which further discussion of constraints on the sense in which such states are symmetric can also be 
found. (In this regard see also Refs. [131 mj.) Here we have demonstrated the symmetry of generic states (not just 
quasiclassical ones) at large volume, in the sense of Eq. (4.12), directly from the exact solutions. 



V. CRITICAL DENSITY AND THE ULTRAVIOLET CUTOFF 



A significant part of the interest in loop quantum cosmology has arisen from the fact that loop quantization seems 
to robustly and generically resolve cosmological singularities; see Ref. [1] [5] for recent overviews. This was noticed 
first in numerical results for semiclassical states [SliSlj subsequently observed in many other models (see e.g. Refs. 



19 



[01 [TD]), and finally proved analytically for all quantum states in the model described in this paper in Ref. "S]. In that 
paper it is shown that the expectation value (and hence spectrum) of the matter density is bounded above by 

^-it-^^^^ (5- la) 

«0.41-pp, (5.1b) 

where pp is the Planck density and the value of the Barbero-Immirzi parameter 7 0.2375 derived from black hole 
thermodynamics has been used |19) . 

We argue here that the existence of a universal upper bound to the density may be traced to the ultraviolet cutoff 
for values of |fc| > |i^/2A| on the eigenfunctions e)^\v). 

Classically, the matter density when the scalar field has value is given by the ratio of the energy in the scalar 
field to the volume, 

In Ref. "3] it is argued that a suitable definition for the corresponding quantum mechanical observable is 

pU = \MI. (5.3) 



where 



Even though the spectrum of this operator is not yet known, an upper bound on the spectrum of A\^ places a bound 
on the spectrum of thence on p\^. Ref. [3j then observed that 

- , (5.5b) 



where \x) is defined through \^) — \/v\4,\x)- Thus, the expectation value of (in the state |^)) may be expressed 
as the ratio of expectation values of the momentum and volume (in the state |x))- Even though the spectrum of 
is not bounded, they go on to show analytically that this ratio is nonetheless bounded above by ^/sJAttj^/ X for all 



states in the domain of the physical observables, leading directly to the bound given by Eq. (5.1 ) on the density. 
Alternately, one might choose to define 



A similar argument then shows 



1 (uj\pI\uj) 

2 {iv\V\i\uj) 



where now j^*) = Fl^jaj). For convenience, we will adopt this latter definition of the density in the sequel. 



A. Heueristic argument 

We offer here a new perspective on the existence of a universal upper bound to the density by arguing that it can 
be seen as a consequence of the linear scaling of the ultraviolet cutoff in the e^*''(i^) with volume, |fc| < |i^/2A|. We 
offer both a new proof of the existence of a critical density in this model in the volume representation, as well as an 
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heuristic argument that has a clear and intuitive interpretation, making it a simple matter to calculate the value of 
the critical density simply from the slope of the scaling of the ultraviolet cutoff. 

In fact, heuristically speaking, using Eq. (5.7) we see that the UV cutoff |fc| < |i^/2A| implies that 
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(5.8a) 



(5.8b) 



(5.8c) 



identical to the rigorous bound on the density - Eq. (5.1) - found in Ref. The linear scaling in the UV cutoff 
on the eigenfunctions thus, in this heuristic way, leads directly to the existence of the universal critical density. In 
particular, the slope of the scaling gives the value of the critical density correctly. 

This "moral" argument is of course not rigorous since and do not commute. There is, however, an interesting 
reason the "moral" argument works: as discussed in Sec. IIIB the scalar momentum and volume operators 



approximately commute, again as a consequence of the ultraviolet cutoff. Thus, also as a consequence of the ultraviolet 
cutoff, the operator is, though its spectrum is not bounded, in effect bounded on subspaces of fixed volume. This 
leads immediately to the upper bound on the density, as in the "moral" argument above. 



The intended meaning of these statements is the following. Because of the ultraviolet cutoff. 
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(5.9a) 
(5.9b) 



Consider the subspace spanned by volume eigenstates with volume less than or equal to some V. While this subspace 
is not strictly invariant under the action of the operator p^, it is approximately so because the off-diagonal matrix 
elements {v\p^\v') are strongly suppressed. The norm of states restricted to this subspace is ||x||y — X]|i/|<v lx(i^)P- 
Then we have 
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(5.10a) 

(5.10b) 

(5.10c) 
(5.10d) 



and we can see that p^ is in effect bounded on subspaces of volume less than a given value. Given Eq. (5.5) or (5.7), 
this helps explain why the heuristic argument above gives the correct value for the critical density. 



B. Proof in the volume representation 



To complete this heueristic argument we offer a new, alternative proof of the existence of a critical density in this 
model in the volume representation, using the definition Eq. (5.6) for the density. (The original proof of Ref. [3] is in 
a different representation of the quantum states and operators and employs the definition Eq. (5.3) for the density, 
though their proof works for either definition.) 

The action of the gravitational constraint in the volume representation, Eq. (|2.9|), may be written 
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(5.11) 
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J2 {'^M'^, 4>) - ^^^{^, 4>r^{^, <t>)} 



is approximately the average of the values of w on either side of the volume v. In this notation, 

{pl)^=h\uomQ\u:{4>)) 

1 fhK " 

2 l2A, 

^u'^uj{v, (j))*uj{u, (j)) 

We wish to show this quantity is bounded above by {Hk / 2Xf' {v\V) ^ ^ in accord with the "moral" argument of Eq 
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(5.8). To proceed, define 
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where v = A\n. Clearly = uj*uj^ so lo and w' have the same norm in the inner product of Eq. (2.11). However, 
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Thus, 
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Since 8 is a positive operator[^we find 
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the sum o: 



Eqs. (5.13) and (5.17) show that the absolute value of t he sum of off-diagonal terms, J2v i^^w*^, is bounded above by 
J the diagonal terms, (i^l?,)!^. Thus, from Eq. (5.13) we see that 
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2A 



(5.18) 



as desired. With the definition Eq. (5.7) for the density, then, the heuristic "moral" argument showing the relation 
between the slope of the scaling of the UV cutoff and the value of the critical density is supported by a direct 
calculation. 



A parallel demonstration in the volume representation using the definition Eq. (5.3) of the density would similarly 
show that the sum of the off-diagonal terms in {Ptt,)x is bounded above by the sum of the diagonal terms. Though 



this can be plausibly argued on the basis of the observations in Sec. IIIB that the off-diagonal elements of in the 
volume representation are bounded above by the diagonal elements, and strongly suppressed for elements connecting 
more than one step off the diagonal - and is of course known to be true because of the proof of Ref. [3 ~ a proof 
entirely in the volume representation at the same level of rigor as that possible for p^ is more difficult because the 
matrix elements of p^ are so much more complicated. The "moral" argument applies in either case. 



^"^ Strictly speaking, since ui' is not generally a solution to the constraint except in regions where w = 0, we should check that G is positive 
on all functions normalizeable in the inner product of Eq. \2.11\ . 
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VI. DISCUSSION 



Working from recent exact results for the eigenfunctions of the dynamical constraint operator in flat, scalar loop 
quantum cosmology, we have demonstrated the presence of a sharp momentum space cutoff in the eigenfunctions 
that sets in at wave numbers \k\ = |j^/2A| that may be understood as an ultraviolet cutoff due to the discreteness 
of spatial volume in loop quantum gravity. Earlier numerical observations showing the onset of a rapid decay in the 
eigenfunctions at small volume at a volume proportional to the eigenvalue ujk are thus confirmed analytically in this 
model. We have argued that the existence of a maximum ("critical") value of the matter density — p'^/'2V\'^ 
that is independent of the state can be viewed as a consequence of the ultraviolet cutoff since the minimum volume 
and maximum momentum scale in the same way. This bound holds for generic quantum states in the theory in the 
domain of the Dirac observables, not only states which are semiclassical at large volume. 

We have offered both an heuristic "moral" argument based on the scaling of the UV cutoff, and a new direct proof 
in the volume representation. While the "moral" argument for the critical density is not rigorous, it is physically and 
intuitively clear, and enables the value of the critical density to be calculated straightforwardly as in Eq. (5.8 1 once 
the slope of the scaling of the cutoff is known. Consistency with the bounds on the matrix elements of the physical 
operators set by the UV cutoff shows the overall coherence of these different points of view. 

It is our hope that this perspective on the origin of the critical density will have some use in the study of more 
complex models. In particular, while the dynamical eigenfunctions have been calculated analytically in this simple 
model, it is probably too much to hope that this will be accomplished in most other, more complicated, models. 
Rigorous proofs of the existence and value of a universal critical density may therefore be difficult to achieve in 
many models beyond sLQC. Nevertheless, in all models it should be possible to study solutions to the gravitational 
constraint numerically. With the recognition from Ref. [3 that the density is bounded by the ratio of the expectation 
value of the momentum to the volume, we have argued here that the existence of a universal critical density may 
be viewed as due to the linear scaling of the ultraviolet momentum space cutoff in the eigenfunctions e\f \v) with 
volume. Therefore, in models in which analytical solutions are not available, numerical evidence for the existence 
of an ultraviolet cutoff in the eigenfunctions may nevertheless be employed to argue robustly for the existence of an 
upper bound to the matter density for generic quantum states in those models, and indeed, its precise value may be 
inferred from the slope of the cutoff scaling. 

The asymptotics enabling the demonstration of the ultraviolet cutoff in the eigenfunctions also enabled us to extract 
analytically the large volume limit of these eigenfunctions based on an analysis of the model's exact solutions. The 
result, consistent with considerable prior work in the field based on physical, analytical and numerical arguments, 
is that the eigenfunctions approach a particular linear combination of the eigenfunctions for the Wheeler-DeWitt 
quantization of the same physical model, with a precise determination of the phase, as well as some understanding 
of the domain of applicability of the approximation. In turn, this allowed us to show that generic quantum states in 
the theory approach symmetric linear combinations of "expanding" and "contracting" Wheeler-DeWitt universes at 
large volume, no matter how non-classical those states may be. 
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